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It is shown through analytical modeling and numerical simulations that external magnetic fields 
can relax the self-trapping thresholds in plasma based accelerators. In addition, the transverse 
location where self-trapping occurs can be selected by adequate choice of the spatial profile of the 
external magnetic field. We also find that magnetic- field assisted self-injection can lead to the 
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be explored using state-of-the-art magnetic fields in current/next generation plasma/laser wakefield 
accelerator experiments. 
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INTRODUCTION 

Plasma based accelerators (PBA) use high intensity 
laser pulses pQ, with intensities above / ~ 10 18 W/cm 2 , 
or highly charged particle bunch drivers [2], with more 
than 10 10 charged particles, to excite ultra-relativistic 
plasma waves. The ideal plasma density to maximize 
charge and energy gain depends on the nature of the 
driver (i.e. lepton, hadron or laser pulse), typically rang- 
ing between no = 10 14 — 10 19 cm -3 . At these plasma 
densities, charged particle bunches can be accelerated by 
plasma wakefields to 1-100 GeV in 1-100 cm [3]. 

The proof-of-principle of PBA is firmly demon- 
strated [4j [5] . Presently, connection with applications [6] 
is an essential step to further improve this technology. 
To this end, fine control over the properties of the ac- 
celerated electrons is required. Several techniques were 
proposed to control self-trapping. Control over charge 
and energy of accelerated bunches can be reached using 
plasma ramps [7 , counter propagating lasers [8] , to ion- 
ization mechanisms [9], and resorting to non-linear opti- 
cal effects such as self-focusing [TO] . 

A novel technique using transverse magnetic fields 
to relax self-injection thresholds has been recently pro- 
posed [11] . The use of external magnetic fields in plasma 
acceleration was first proposed to extend the accelera- 
tion distances in plasma accelerators in the surfatron 
model [12]. The role of external magnetic fields in PBA 
was also explored in [13] , and the use of longitudinal mag- 
netic fields to enhance the self-injected charge in laser 
wakefield acceleration was investigated in [14 . 

This paper presents a detailed derivation of the self- 
trapping threshold condition in the presence of external 
fields. Using the particle-in-cell (PIC) code Osiris [15], 
it is shown that magnetic injection can be used to gener- 
ate single or multiple off- axis self-injected bunches with 
well defined radial injection positions. Using the post- 



processing radiation code JRad [16] it is demonstrated 
that these electrons may emit betatron radiation at well 
defined frequencies close to the undulator regime. This 
paper is structured as follows. In Sec. | describes an 
analytical trapping condition in the presence of external 
fields. In Sec. [] 3D PIC simulation results are employed 
to analyze the relevant physical mechanisms of magnetic 
field assisted self-injection. The use of different B-field 
geometries to control transverse features of magnetically 
injected electrons is described in Sec. [] Section [] shows 
that magnetically assisted injection can lead to the emis- 
sion of clearly defined betatron radiation harmonics for 
the first time in PBAs. Finally, conclusions are stated in 
Sec. D 



TRAPPING CONDITIONS IN THE PRESENCE 
OF EXTERNAL FIELDS 

The dynamics of the electrons in the fields created 
by an intense laser in the blowout regime can be de- 
scribed using Hamiltonian dynamics [9]. A general trap- 
ping condition in the presence of external fields can be 
found by examining the evolution of the Hamiltonian of 
plasma electrons in the co-moving frame, {x = x : y = 
U->i — v (f>t — z,s = z), given by H = H — v^Pn, where 
(x, y) are the transverse coordinates, z the distance, 
the wake phase velocity (determined by the driver 
group velocity), P\\ the longitudinal canonical momen- 
tum, H = ^/ra^c 4 + (P + eA) 2 — ecj) the Hamiltonian of 
a charged particle in the presence of electric and magnetic 
fields, m e and e the electron mass and charge, c the speed 
of light, A and <p the vector and scalar potentials, and 
P = p — eA, where P and p are the canonical and linear 
momentum respectively. Normalized units will be used 
henceforth unless explicitly stated. Mass and charges are 
normalized to m e and e, respectively, velocity v to c, time 
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to uj p = Y // 4 / 7rnoe 2 /m e , momentum to ra e c, and density 
to the background plasma density no- Vector and scalar 
potentials are normalized to e/m e c 2 and e/m e c respec- 
tively. Magnetic fields (B) are normalized to uj c /oj p where 
uj c = eB/m e is the cyclotron frequency. 

In order to derive a trapping threshold condition in the 
presence of external fields, we consider first the expres- 
sion for the temporal evolution of H: 
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Integration over the particle trajectory yields: 
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where the subscripts V and re f er to the initial and 
final (trapped) electron positions. The integration is per- 
formed along the electron trajectory, and dt = d£/(i>0 — 
v z ). Combining Eq. ([!]) with Eq. Q gives: 
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Using the definition for H, and considering that initially 
electrons are at rest (i.e. = 0): 

Uf-Ui =7/-W||-l-(0/ -VjAf^-fa -v^w) , 

(4) 

with <j> = ^ + ^ ext , A\\ = Af + Af\ and where the 
superscripts 'pi' and 'ext' refer to the plasma and external 
fields respectively, and 7 = (1 — v 2 ) -1 / 2 is the relativistic 
factor. Defining the wake potential i/j = <fi — v<j,A\\ , = 
i/jf — ^i, and assuming that for trapping the longitudinal 
velocity of the electron must reach the velocity of the 
wake i.e. v z = Eq. Q readily becomes: 
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where 7^ = (1 — v 2 ^)~ x l 2 is the gamma factor of the 
plasma wave. Using Eq. ([3| to express Hf — Hi leads 
to the trapping condition [11] : 
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Equation (|6| is a general trapping condition in the pres- 
ence of external fields, and valid beyond the range of va- 
lidity of the quasi-static approximation [17]. Analytical 
solutions to Eq. (J6|, however, are not yet known because 
the calculation of A^ pl and J d^Hdt; requires accurate 
prediction of the particle trajectories and field structures 
at the back of the bubble, where the applicability of the 
standard analytical models [18] is limited. 

To retrieve a general trapping threshold in the absence 
of external fields, and in the conditions where the quasi- 
static approximation is valid, it should be considered 



A ^ext = and j ^ftdf = in Eq. For an ultra- 
relativistic plasma wave with 7^ — > 00 trapping occurs 
when 1 + Aip pl = 0, or, equivalently, A^ pl = — 1. Gen- 
erally, this condition can only be met at the back of the 
plasma wave, in regions of maximum accelerating fields, 
where is minimum and approaches — — 1 [9]. In 
the presence of static external fields the trapping condi- 
tion becomes 1 + A^ pl = — Atp ext . The trapping thresh- 
olds are relaxed because they can be met when Aip pl is 
larger than —1 provided that A^ ext > 0. In other words, 
trapping may occur for lower values of peak accelerat- 
ing gradients. Moreover, trapping may be suppressed if 

A ^ext < q 

Trapping can also be relaxed (or suppressed) when the 
external fields vary spatially in z because of the contri- 
bution of finite / d^Hd£ ^ to Eq. If the pro- 
file of the external fields profile leads to J d^Hdt; > 
(J d^Hd^ < 0) along the electron trajectory then trap- 
ping is facilitated (suppressed) [10] [TT] . Physically, the 
fact that J d^Hdt; > is typically associated with the re- 
duction of the wake phase velocity through the accordion 
effect, thus facilitating self-injection [TO] [TT] [T9]. 



MAGNETICALLY CONTROLLED 
SELF-INJECTION IN LWFAS AND PWFAS 

To investigate controlled self-trapping in the presence 
of external static magnetic fields we present in this sec- 
tion 3D particle-in-cell simulations of laser (LWFA) and 
plasma (PWFA) wakefield accelerators using the particle- 
in-cell code Osiris. Figure [T] illustrates the evolution, 
and highlights the key mechanisms of injection assisted 
by external B-fields in the LWFA. The simulation win- 
dow moves at the speed of light, with dimensions of 
24 x 24 x 12 (c/uj p ) 3 , divided into 480 x 480 x 1200 
cells with 1x1x2 electrons per cell in the (x, z) 
directions respectively. The plasma ions are immobile. 
A linearly polarized laser pulse with central frequency 
uoq/uo p = 20 was used, with a peak vector potential of 
ao = 3, a duration cj p tfwhm = 2-^/oq, and a trans- 
verse spot size matched to the pulse duration such that 
Wo = ctfwhm [3]. The plasma density is of the form 
n = no(z) (l + Anr 2 ) for r < y/W c/uo p and n = for 
r > y/l0c/uj p with An = An c = 4/Wq (i.e. the nor- 
malized matching condition given by An c = 4/(7rr e WQ ), 
where r e = e 2 /m e c 2 is the classical electron radius) be- 
ing the linear guiding condition in the normalized units, 
and where uq(z) is a linear function of z which increases 
from no = to no = 1 for 50 c/uj p to ensure a smooth 
vacuum-plasma transition. The channel guides the front 
of the laser thereby minimizing the evolution of the bub- 
ble. A static external B-field pointing in the positive 
y-direction was used. At the point where the plasma den- 
sity reaches its maximum value, the external field rises 
with = uj c /uj p = 0.6sin 2 [7rz/(2L ramp ) + $1], with 
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Z/ amp = 10c/ 'u p . It is constant and equal to = 0.6 
for L flat = 40 c/ujp and then drops back to zero with 
B* xt = 0.6sin 2 [7r^/(2L ramp ) + $ 2 ]- Moreover, $i and 
<l>2 are phases chosen to guarantee the continuity of the 
external B-field profile. Qualitatively, the longitudinal 
profile of the magnetic field thus resembles that of [20 . 
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FIG. 1: 3D osiris simulation results illustrating the mag- 
netic self-injection mechanisms, (a), (c), and (e) show the 
electron plasma density in gray, the self-trapped particles in 
blue, and the laser pulse envelope in red colors at t = 110 /uj p , 
t = 126/cjp and t = 159/cj p . (b), (d), and (f) show the cor- 
responding p|| — £ phase-space. The magnetic field leads to 
off-axis self-injection in a narrow angular region. The inset 
in (f) represents the transverse momentum phase- space of the 
self-injected electron bunch residing within the bubble. The 
B-field profile is schematically represented on the top of the 
figure. The laser driver moves from left to right as indicated 
by the arrow. 

Self-injection is absent from the regions where the B- 
field rises. In these regions, electrons traveling backwards 
(v z < 0) near the back of the bubble feel an increas- 
ing v x B ext force that rotates electrons anti-clockwise 
thereby locally decreasing (increasing) the blowout ra- 
dius for x > (x < 0). Then, as the B-field rises, the 
local wake phase velocity at the back of the bubble in- 
creases (decreases) for x > (x < 0). For x > 0, 
is superluminal, J d^Hd^ < 0, and self-injection can not 
occur. For x < trapping is precluded because electrons 
reach the axis in regions where the plasma focusing and 
accelerating fields are unable to focus and trap electrons. 
Thus, although for x < J d^Hd^ > 0, we have that 
/ d^ftdf + A^ ext < 0. 

Self-trapping occurs in the uniform regions of the ex- 
ternal magnetic field where x > 0. For x > 0, elec- 
trons rotating ant i- clockwise reach the axis in regions of 
maximum focusing and accelerating fields with larger py 
and can be trapped. For x < 0, electrons reach regions 
of the axis (where focusing and accelerating fields are 
lower) with lower py, and are are lost to the surround- 
ing plasma. A threshold B-field for injection may be 
retrieved in the limit where 7^ — oo. Neglecting the 
plasma fields (A^ pl = 0), and noting that the exter- 
nal longitudinal vector potential 74jj xt = —B y x is consis- 
tent with the considered magnetic field, leads to a sim- 



plified trapping condition A?/; ext = —B y Ax — 1, where 
Ax = Xf — Xi ~ — 7*b, where is the blowout radius and 
where it was considered that the initial (final) trapped 
electron radial position is x = (x = 0). It shows that 
injection is facilitated in the region where Ax < is min- 
imum. As a consequence, injection occurs off-axis (for 
x > 0), and in a well defined azimuthal region defined by 
— ByTb sin 6 = 1, where is the angle between the plane 
of the electron trajectory with the B field [11]. Note, 
however, that this trapping threshold condition overes- 
timates the threshold B-field for self-injection because it 
neglects the plasma fields. 

There is an upper B y value, given by uo c /uo p < 1, be- 
yond which injection may be suppressed in regions where 
the B field is flat. The later condition ensures that the 
plasma wakefields are nearly undisturbed by the exter- 
nal fields. Simulations then showed that when u c /ujp ^> 1 
there is a suppression of the wakefields that prevents in- 
jection. Hence trapping can be relaxed in the regions of 
uniform B-fields provided that 1/r^ < B y < 1 or, equiv- 
alently 170/r 6 [10/im] < B[T] < 32 v /n [10 16 cm- 3 ]. 

The above-mentioned upper B-field limit for injection 
is absent from the downramp regions, where a stronger 
self-injection burst occurs for x > 0. Injection occurs 
within the same angular and radial region as in the uni- 
form B-field section (Fig. [l^). For x > 0, when the B- 
field lowers increases, lowers and J d^Hdt; > 0, 
facilitating injection. For x < 0, v$ > 1, and trapping is 
suppressed. The resulting phase space at t — 126c/ uj p is 
shown in Fig. [T]i. 

After the magnetized plasma region, the magnetically 
injected electron bunch is clearly detached from the back 
of the bubble, leading to the generation of a quasi- mono- 
energetic electron bunch. The magnetic injected electron 
bunch right after the B-field is shown at t = 159c/ 'uj p 
in Fig. and the corresponding phase-space in Fig. [if. 
The inset of Fig. [if shows the transverse phase space of 
the magnetically injected electron bunch residing within 
the blowout region. The asymmetrical distribution re- 
sults from the fact that the injection process occurs off- 
axis. At this location, the emittance of the beam is on 
the order of Itt mm mrad in both transverse directions. 
Although comparison of beam emittance with a similar 
scenario without the B-field is not meaninfull because 
without the field the amount of self-injected charge is 
much smaller. However, the measured beam emittance 
is at the same values or lower than typical emittances of 
LWFAs. 

External magnetic fields also relax the self-trapping 
thresholds in the PWFA. Figure [2] shows results from a 
3D simulation of a magnetized PWFA. A 30 GeV elec- 
tron bunch was considered with density profile given 
by n b = n 60 exp (-x 2 L /(2a 2 1 _ )) exp (-£ 2 /( 2cr f )) > with 

<j± = 0.17 c/cjp, <j£ = 1.95 c/cjp, and n^/no = 19. These 
parameters ensure that is similar to the magnetized 
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LWFA investigated above. The simulation window di- 
mensions are 24 x 24 x 24 (c/cj p ) 3 , and it is divided in 
480 x 480 x 640 cells with 1x1x2 electrons per cell in 
the (x, y, z) directions respectively. The magnetic field 
profile is similar to the LWFA case. 
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FIG. 2: 3D osiris simulation results illustrating magnetic self- 
injection in the plasma wakefield accelerator using (a) B eyit — 
0, (b) £ ext = 0.2, and (c) B ext = 0.6 at t = 275/u p . In this 
scenario the B-field injection mechanism is dominated by the 
plasma bubble dynamics in the B-neld down-ramp. It is clear 
that the trapped charge increases with the amplitude of the 
external field. The driver moves from left to right as indicated 
by the arrow in (c). 



The inset of Fig. ^ also shows that the energy spreads 
of the magnetically injected electrons are on the order 
of 100%. Due to the short duration of the self-injected 
bunch in comparison to the plasma wavelength, which 
guarantees uniform acceleration throughout the entire 
bunch length, the relative energy spread would decrease 
as the beam accelerates. Moreover, the energy spread 
would further narrow down near the dephasing length 
due to the bunch phase-space rotation [22 . 

For these parameters the threshold magnetic field for 
self-injection is By Xt > 0.2. To connect these sim- 
ulations with actual experimental conditions, we take 
no = 10 1T cm -3 for which the electron beam and plasma 
parameters match those available at SLAC [5] with a± = 
50.4 /im, a z = 84 fim and a total number of 3 x 10 10 elec- 
trons. For these parameters, = 0.2 corresponds to 
20 T. These magnetic fields could be produced with state- 
of-the-art magnetic field generation techniques [2Q1 [25] . 
By tuning further the plasma parameters controlled in- 
jection with magnetic fields as low as 5 T can also be 
achieved (cf. Sec.[|. 



The accelerating structures are similar for the LWFA 
and PWFA parameters described above since the blowout 
radius are similar for both cases. However, as shown by 
Eq. (J6|, self-injection thresholds are harder to meet in the 
PWFA than in the LWFA because 7^ WFA ~ 60000 > 
^lwfa 20. In contrast to the LWFA scenario, injec- 
tion is then absent in the PWFA in the uniform regions 
of the B-field, where the larger p z at the back of the bub- 
ble for x > 0, associated with the additional electrons 
v x B anti-clockwise rotation, is still below that required 
for injection. Magnetic self-injection occurs only in the 
B-field downramp (Fig.j2j>c), where the injection mecha- 
nism is similar to that ascribed to the LWFA. In general, 
stronger self-injection bursts occur in the B-field down- 
ramp for both LWFA and PWFA. 

The amount of self-injected charge can be tuned by 
changing the B-field amplitude. The inset of Fig. (2^ 
shows the spectra of the self-injected charge in the first 
plasma bucket using B ext = 0.6 (red curve), B ext = 0.2 
(green curve), and £? ext = 0.0 (blue curve). The amount 
of trapped charge is negligible in the un-magnetized sce- 
nario, and it is roughly 8 times larger for £? ext = 0.6 
than for B ext =0.2 (notice that the plot is logarithmic in 
the vertical y-direction). These results show that higher 
B-field amplitudes increase the total amount of injected 
charge. 

Because of beam- loading [21] , higher amounts of self- 
injected charge lead to lower accelerating gradients. Con- 
sequently, the maximum energy that can be achieved is 
lower for self-injected bunches with higher charges. This 
is consistent with the inset of Fig. ^jp which shows that 
self-injected bunches with lower charges reach higher en- 
ergies. 



SIMULTANEOUS GENERATION OF MULTIPLE 
SELF-INJECTED ELECTRON BUNCHES 

The transverse location where self-trapping is relaxed 
can be selected by adequate choice of the profile of the 
external magnetic field. As an example, Fig. [3] shows the 
results from a 2D slab geometry simulation using a mag- 
netic field which reverses sign at x = (this is equivalent 
to an azimuthal B-field profile in cylindrical symmetry). 
In this case, the magnetic field points outside (inside) 
the simulation plane for x > (x < 0). The 2D sim- 
ulations use a simulation box that moves at c with di- 
mensions 12 x 32 (c/ujp) 2 , and is divided into 640 x 3000 
cells with 3x3 electrons per cell in the (x, £) directions 
respectively. The laser pulse and plasma channel param- 
eters are similar to those of the 3D LWFA simulation 
(cf. Fig. [I]). The amplitude of the external B-field is 
£* xt = uj c /uo p = 0.6sin 2 [7rz/(2L ramp ) + $i]x/\x\, with 
L ramp = 10c/cjp, it is constant and equal to By Xt = 0.6 
for L flat = 50 c/ujp and drops back to zero with By Xt = 
0.6sin 2 [7rz/(2L ramp ) + &2]x/\x\, where the choice of $ x 
and $2 ensures the continuity of the B-field longitudinal 
profile. 

Figure^ shows the magnetically injected electrons in 
the regions where the B-field is uniform. Two off-axis in- 
jection bursts occur at well defined transverse positions 
in the fiat B-field regions. The two bunches are then in- 
jected symmetrically close to x = 0. An additional and 
stronger self-injection burst occurs at the B-field down- 
ramp (Fig. [3}d). After the magnetized plasma region, the 
two self-injected electron bunches continuously acceler- 
ate in the wakefield (Fig. [3J3). Note that Fig. [3J3 refers 
to the early propagation of the electron bunch, much 
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shorter than the dephasing length. Similarly, the propa- 
gation distance is much smaller than the betatron period 
of oscillation. The physical mechanisms under which self- 
injection occurs in the present configuration are identical 
to those presented in Sections [] and [J 

Interestingly, Fig. [3] reveals that injection occurs in a 
highly spatially localized region. Off axis injection from 
well defined radial and azimuthal regions was observed in 
Sec. [] in 3D simulations. Generally, however, this effect 
is more noticeable in 2D slab geometry simulations than 
in 3D. These results also suggest that ring like electron 
bunches could be obtained in 3D. This could be advan- 
tageous for radiation generation purposes because bunch 
particles would perform betatron oscillations with similar 
amplitudes. 
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FIG. 3: 2D Osiris simulation of a magnetized LWFA using an 
azimuthal-like magnetic field. In the upper (lower) simulation 
mid-plane the B-field is directed outside (inside) the simula- 
tion plane, (a) is taken in the B-field flat region at t = 79 /uj p , 
(b) in the B-field down-ramp at t = 92/cj p , and (c) after the 
B-field at t — 102 /uo p . The B-field profile is schematically 
represented on the top of the figure. This configuration leads 
to identical magnetically injected electron bunches in the up- 
per and lower simulation mid-plane. The laser driver moves 
from left to right as indicated by the arrow. The dark points 
represent self-injected particles. 



EMISSION OF BETATRON RADIATION AT 
WELL DEFINED FREQUENCIES 

Typical synchrotron radiation experiments in plasma 
accelerators reveal that radiation emission occurs in the 
wiggler regime. The wiggler regime enables emission of 
x-rays with broad spectra [23]. This contrasts with the 
undulator regime, where radiation is emitted at well de- 
fined harmonics. Although not yet attained experimen- 
tally, the undulator regime provides ideal conditions for 
radiation amplification, being critical for the realization 
of a ion-channel plasma based laser [24]. This section il- 
lustrates how could magnetically injected electrons emit 
betatron radiation at well defined frequencies, closer to 
the undulator regime. 

The PWFA beam and plasma simulation parameters, 
presented in Sec. [J were tuned in order to lower the re- 
quired magnetic field for injection, such that it could be 
more easily reached experimentally, and in order to lower 



the amplitude of the betatron oscillations in comparison 
to the plasma skin depth, such that distinguishable be- 
tatron radiation harmonics could be emitted. System- 
atic 3D parameter scans then showed that the threshold 
magnetic field for injection is 5.5 T at no = 10 15 cm -3 . 

mm, and 



At this plasma density, L 

L ramp = 1Q ^ 



flat 



40c/uj p 



6.S 



1.68 mm, and the maximum B-field 
amplitude is = 0.55 uj c /uo p . These parameters are 
within current technological reach [20j [25] . Simulations 
used a simulation box with 12 x 12 x 16 (c/uj p ) 3 , divided 
into 480 x 480 x 640 cells with 2x2x1 particles per cell 
for the electron beam and background plasma. 
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FIG. 4: (a) B-injected electron bunch trajectories in the 
plasma wakefield accelerator. The trajectories gentle rise in 
the X2 direction is due to the subtle transverse velocity shift 
imprinted on the particle beam driver in the B-field region. 
The inset represents the distribution of the radiation strength 
parameter ap. A significant fraction (83%) of the trajectories 
are characterized by ap <C 1. (b)-(c) Corresponding radiation 
profile lineouts passing through the center of a virtual detector 
located 5100c/ oj p from the end of the plasma. The white 
dashed lines corresponds to the predictions of Eq. |7| using 
7 = 400 and rp = 0.06. The location where the spectra were 
taken in the detector plane is also shown. 



Self-trapping occurs off- axis at the B-field downramp. 
Injection is localized radially and azimuthally, enabling 
the bunch to perform synchronized betatron oscillations 
(Fig. [4^l). The transverse x-axis shift of the electron tra- 
jectories result from the electron beam driver deflection 
when traversing the magnetized plasma region. The de- 
flection angle is small and could be corrected by adding 
additional magnetized plasma regions with alternating 
B-fields along the propagation direction [TT] . 

The small blowout radius (r& ~ l.bc/uj p ) ensures 
that the betatron amplitudes of oscillation rp are much 
smaller than the plasma skin-depth ((rp) c± 0.06c/uj p ). 
The corresponding radiation strength parameter (ap = 
^Kprpkp) distribution, where Kp = l/^/2j is the nor- 
malized betatron frequency, and k p the plasma wavenum- 
ber, is shown in the inset of Fig. [4^,. It shows that a 
significant portion (83%) of the electrons radiate with 
ap < 1, an indication that single harmonics could be 
distinguishable in the emitted radiation spectrum. To 
retrieve the radiation spectrum, a random sample of the 
self-injected electrons was post-processed using the radi- 
ation code JRad [16 . Figure j4j> c shows the radiation 
spectrum in the transverse central lines of a virtual de- 
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tector placed at a distance 5100 c/u p from the exit of the 
plasma. The detector lies on the x — y plane. 

The asymmetries in the x direction depicted in Fig. 
result from the tilt of the magnetically injected electron 
trajectories. Figure ^p-c reveal that radiation is emitted 
at well defined frequencies, which are particularly clear 
at larger angles, i.e. for larger \x\. The width of each 
harmonic present in Fig. is larger than that expected 
in an idealized scenario, where radiation would be purely 
emitted in the undulator regime. This widening is due 
to the spread on the ap distribution (through 7 and rp 
spreads) and also because some electrons radiate with 
strength parameters which are larger than unity ap > 1. 

For an electron bunch with constant relativistic factor 
7, and constant rp in a pure ion-channel, the frequency 
of the betatron radiation harmonics emitted in the un- 
dulator regime are given by [26] : 

= 2nj 2 Kp 

u p (^l + a2/2^cos6> + 27 2 (l-cos(9)' 

where n corresponds to the n th emitted harmonic, and 6 
to the angle between the velocity vector of the electron 
and the point in the detector. Radiation collected on- axis 
only exhibits odd-harmonics. To compare the predictions 
of Eq. ([7]) with simulation results we computed the parti- 
cles trajectories average (7) = 400 and (rp) = 0.06. This 
yields ap ~ 0.7, consistent with the inset of Fig. [4j The 
analytical prediction Eq. ^ are shown by the dashed 
lines in Fig. j4j>c. Eq. ^ is in good agreement with the 
simulation results specially for larger values of \x\. Dis- 
crepancies are due to the fact that the beam trajectories 
are tilted, and that ap, rp, and 7 vary in time and for 
each electron. 



CONCLUSIONS 

In conclusion, we explored further a recent controlled 
injection technique that uses transverse, static mag- 
netic fields to tailor transverse properties of self-injection. 
This scheme leads to off- axis self-injection in well de- 
fined radial and azimuthal regions. A configuration con- 
sisting of a section of transversely uniform magnetized 
plasma yielding off- axis self-injection was investigated. 
It was shown that simultaneous self-injection of electron 
bunches could be achieved by using transversely non- 
uniform fields. This work also suggests that a series of 
magnetized regions could be used to produce a temporal 
sequence of electron bunches. Moreover, multiple spa- 
tially separated electrons could be produced simultane- 
ously with transversely non-uniform B-fields. We showed 
that this technique could be used to produce electron 
bunches capable to emit betatron radiation at well de- 
fined frequencies with current technology. 
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